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Abstract. A family of orthonormal polynomials on the unit ball B"* of ] 
with respect to the inner product 

(/, g)= [ A[{1 - \\xf)f (x)]A[(l - \\x\\)g(x)]dx, 

where A is the Laplace operator, is constructed explicitly. 



1. Introduction 

In a recent study |lj on the numerical solution of the nonlinear Poisson equation 
—Am ~ f{',u) on the unit disk with zero boundary conditions, Atkinson PP asked 
the question of finding an explicit orthogonal basis for the inner product defined by 

(/, 5)a := - / A[(l -x^- y^)f{x, y)] A[(l - - y')g{x, y)]dxdy 

on the unit disk of the Euclidean plane, where A is the usual Laplace operator. 
The purpose of this note is to provide an answer for this question. 

We shall consider more generally the analogous inner product on the unit ball 
B"^ in Mf^. We call orthogonal polynomials with respect to such an inner product 
Sobolev orthogonal polynomials. In the theory of orthogonal polynomials of one 
variable, the name Sobolev is associated with polynomials that are orthogonal with 
respect to an inner product defined using both functions and their derivatives; see, 
for example, 0] and the references therein. As far as we know, Sobolev orthogonal 
polynomials have not been studied in the case of several variables. 

Our main result, given in Section 2, is a family of orthonormal polynomials with 
respect to (•, •)a on B"^ that are constructed using spherical harmonics and Jacobi 
polynomials in Section 2. For d = 1, orthogonal polynomials with respect to this 
inner product has been studied recently in The explicit formula can be used 
to study further properties of the orthogonal basis. In particular, it turns out that 
the orthogonal expansion of a function / in this basis can be computed without 
involving the derivatives of /. This will be discussed in Section 3. 

2. Sobolev Orthogonal Polynomials 

For X G M'*, let ||a;l| denote the usual Euclidean norm of x. The unit ball in R*^ 
is B'' := {x : \\x\\ < 1}. Its surface is S'^'^ := {x : \\x\\ = 1}. The volume of B'^ 
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and the surface area of S"^ ^ are denoted by vol(i?'^) and Wd-i, respectively, 

vol(B'') = iva-i/d and uja-i = 271^^/^ /r{d/2). 

Let 11'' = M[xi, . . . , Xd] be the ring of polynomials in d variables and let denote 
the subspace of polynomials of total degree at most n. We consider the inner 
product defined on the polynomial space by 

if^9)A / A[(l - ||xf )/(x)]A[(l - \\xr)gix)]dx. 



4d2 vol(B'') Jb 

The constants are chosen so that (1,1) a — 1- As pointed out in [T], the inner 
product is well defined and positive definite on 11''. Let V^(A) denote the space of 
orthogonal polynomials with respect to (/,<?) a- It follows from the general theory 
of orthogonal polynomials in several variables (jSI) that the dimension of V^(A) is 
("+^7^). A polynomial P is in V^(A) if it is orthogonal to all polynomials of lower 

degrees with respect to (/, .g)A- If {-Pa} is a basis of V^(A) and {Pa,Pp)A = 
whenever a 7^ /?, it is called a mutually orthogonal basis. If, in addition, Pa is 
normalized so that {Pa, Pa) A — 1, the basis is called orthonormal. Our objective 
in this section is to find an explicit orthonormal basis for V^(A). 

The presence of the Laplace operator suggests that we make use of harmonic 
polynomials, which are homogeneous polynomials that satisfy the equation AP — 0. 
Let denote the space of harmonic polynomials of degree n. It is well known 
that 

^^'"^"^i d-1 )-[ d-1 

The restriction of y e Hf^ on S"^^^ are called spherical harmonics. They are 
orthogonal on S"^^^. We will use the spherical polar coordinates x — rx' for x G M'', 
r > 0, and x' £ S"^-^ . For Y G we use the notation Y[x) to denote the harmonic 
polynomials and use Y{x') to denote the spherical harmonics. This agrees with 
X = rx' since y is a homogeneous polynomial, Y{x) = r'^Y{x'). Throughout this 
paper, we use the notation {Y^ : 1 < < (T„} to denote an orthonormal basis for 
Tifj, that is, 

(2.1) f y;(x')yr(a;')dc^(x') = 5^.Jn,m, 1 < M, ^ < fT„, 

where dcu stands for the surface measure on S'^~^. In terms of the spherical polar 
coordinates, x = rx', r > and x' S S'''"^, the Laplace operator can be written as 

. d-1 d 1 ^ 

2.2 ^ = T^ + ^ + -^0, 

or^ r or 

where Aq is the spherical Laplacian on 5''^^. It is well-known that 

(2.3) Aor(.T') = -n{n + d- 2)Y{x'), Y G H^, x' G S'^-^. 

The spherical harmonics have been used to construct orthogonal polynomials on 
the unit ball. For later use, let us mention an orthogonal basis with respect to the 
inner product 

(/,.9)p / f{x)g{x)W^{x)dx, W^{x) = (1 - \\xfY, 
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where /i > — 1 and is the normahzation constant of W^. Let V,'f(PF^) denote 
the space of orthogonal polynomials of degree n. A mutually orthogonal basis for 
V„(W^) is given by (PI) 

(2.4) PlAW,:x) = P,('^'"^'^+^\2||xf - l)Yr^^{x), 0<j< n/2, 

where pj-"'^-' denotes the Jacobi polynomial of degree j, which is orthogonal with 
respect to (1 — x)°'{l + x)^ on [—1, 1], and {Y^~^^ ■ ^ < j ^ <^n-2j} is a basis for 

n-/d 

''■n-2j- 

In view of 12.41) we will look for a basis with respect to (/, .g)A in the form of 

(2.5) Ql,{x)^q,{2\\xf-l)Yr'Hx), 0<2j<n, Y^^^ eUt^^^, 
where qj is a polynomial of degree j in one variable. 

Lemma 2.1. Let Q"^ he defined as above. Then 

A [(1 - \\xr)Ql,ix)] =HJpq,) (2r2 - l)Yr''ix), 
where (3 — n — 2j + and 

{JpQj)is) = (1 - s')q^{s) + (/3 - 1 - (/? + 3)s)q'^{s) -{(3 + l)q,{s). 
Proof. Using spherical-polar coordinates, H2.2|l and l|2.3|l show that 
A [(1 - \\xr)Ql,{x)] = A [(1 - r')q,i2r' l)r"-2.yn-2.(^')] 

= 4r"-2, [4^2(^ _ r^)q';{2r^ - 1) + 2((/3 + 1) - (/3 + '3y)q'j{2r^ - 1) 
-{P+l)q,{2r^-l)]Yr''{x'). 
Setting s 2r^ — 1 gives the stated result. □ 
Lemma 2.2. Let e Ilfc :— 11^ be orthogonal with respect to the inner product 

if,9)p-^ fypfKs)iJpg)is)il + sfds, P>-1. 

Then the polynomials Qj\, in (12. 5|) with qj = Pj"'^^^ where Pk = k + {d — 2)/2, 
form a mutually orthogonal basis for V,^(A). 

Proof. It is easy to see that {f,g)i3 is indeed a positive definite inner product on 
the space of polynomials of one variables, so that the orthogonal polynomials with 
respect to {f,g)i3 exist (see Lemma |01 below). Using the formula 

f{x)dx^ [ r'^-^ [ f{rx')duj{x')dr, 



the definition of and H2.1|) shows immediately that 

X ^^V'^+'("-'^'^-^42(Jft._.,g,)(2r^ - l)(J/3„,_,,,<Z/)(2r2 - l)dr. 

In the nonzero case we have (3n-2j — Pn'-2j'. Thus, a change of variable r ^-^ 
^(1 + s)/2 shows that 

(2-6) (<3j,i/' Qf.y)A = S^ySn-2j,7i' -2j' -j^p:;^{qj , q]')p^-2j : 

which proves the stated result. □ 
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We note that ^' should be understood as one member (of degree j) in the 
orthogonal family {^f""^^}- 

Lemma 2.3. The polynomials Pj defined by 

pUs) = 1, p';{s) = il-s)pj'jpis), j>l, 
are orthogonal with respect to the inner product {f,g)fj. 
Proof. We need the following property of the Jacobi polynomials fHl, p. 71], 



(2.7) 



(l-s)F|!f (s) = 



23+13+1 



{j + l)P^'+P{s)-jP)''^'{s) 



(l,/3)/ 



The Jacobi polynomial pjl'f satisfies a differential equation 

(1 - s^)y" + (-1 + /3 - (3 + P)s)y' + [j - l)(j + /3 + l)y = 0. 
Using these two facts, we easily deduce that 

'-^-^J, [(1 - s)Pl':?\s)\ . (, + DJ.P^'jpis) - jJ,P^'J\s) 
= (j + 1) [i-U m + /3 + 1) - (/3 + (s) 

-j{-j{j+f3 + 2)-{P + l))P^''^\s) 
= -jU + 1) [{j + l3)P^'-?\s) - {j +13+ l)Pj'^^'\s) . 
We need yet another formula of Jacobi polynomials 2, p. 782, (22.7.8)], 

(2.8) (2j +(3+ l)Pj"^^'\s) ^{j+P+ l)Pj'^^'\s) - (j + P)Pj'J\s) 
which implies immediately that 

(2.9) J, [(1 - s)P^'J\s)] = 2j(j + l)Pj°-''\s). 
Hence, for > 1, we conclude that 



(2.10) (pf,pf,)/3 = / j0 ii-s)p^'j^'is) J, (i-s)4rAs) {i + ^r)d. 



(2,/3)/ 



3(2,/3), 



m + !)/(/ + 1) / Pf^\s)Pf^\s){\ + sfds = 



whenever j ^ j' . Furthermore, for j > 1, we have 

{P^„P^)0 = -2j(j + l)(/3 + 1) Pj"^^'\s){l + sfds = 
since iJpp^^)is) = iJpl) = -(/3 + 1). 



□ 



As a consequence of the above lemmas, we have found a mutually orthogonal 
basis with respect to (•, •)a- 

Theorem 2.4. A mutually orthogonal basis for V,f (A) is given by 

(2-11) ^„ , , ,^ „ „2,„(2,n-2j + ii^),.„ o., , . _ . n 



o;^,(x) = (i-iix||^)p;_r 
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where {Y^ ■.\<v < cr„_2j} is an orthonormal basis ofTi,n_2j- Furth 
(2.12) {Qo,u^Qo,u)a = — (<3",y,Q"^)A = "' '-■^'^ -* 



d{n + d/2) 



Proof. The fact that Q'l^ e V^(A) foUows from Lemma EH and Lemma lO To 
compute the norm of Qq ^ we use the fact that 

(2.13) A[(l - \\xr)Yr'Hx)] - -2d>^;(x) -4(x,V)r," - -2(d + 2n)i;"(a;) 
by Euler's formula on homogeneous polynomials, which shows that 



Furthermore, using the equation H2.5|l and H2.10|l . we have 

d{Pj+2j + l) d{n + d/2) 

where we have used the well known formula for the norm of the Jacobi polynomial 
(see, for example, [HI p. 68]). □ 

The explicit formula of the basis H2.11|l leads to the following interesting result, 
which relates V,'J(A) to orthogonal polynomials with respect to W2{x) = (1 — ||a;||)^. 

Corollary 2.5. For n > 1, 

V^l{A)^nt(Bil^\\xf)Vt,iW2). 
Proof. Using the basis (|2.4|l for 1(14^2), it follows that wc actually have 
(2.14) Q]Jx) = (1 - \\xr)P;^,'jW2;x) 

for j > 1, from which the stated result follows. □ 
In the case of = 2, an orthonormal basis for the space is given by 

r" (x, y) = -y/ir" cos ne and Y^' {x,y) = y^r" sin nO 

in polar coordinates x — r cos 9, y — r sin 6. Hence, a mutually orthogonal basis for 
V,^(A) is given by 

Ql,{x,y)=Y,^{x,y), Ql^ix , y) ^ Y^ {x , y) , 

Ql,{x,y) ^{l~x^^ y^)p\^_:r^'\2x^ + 2y^ l)Yr'\x,y), 1 < j < f 

Ql2{x,y) ^{l-x^- y')P^','r''\^^' + V - l)Yr''{x,y), 1<J<^, 

which becomes an orthonormal basis upon dividing by the square root of the norm 
given by 12.12|l . Without normalization, this gives 

Vf{A) = span{2;,2/}, V|(A) = span{x^ ~ y^ ,xy,l ~ x^ ~ y^}, 

V|(A) = span{x3 - Sxy^, Sy^ - x^y^ x{l - x^ - y^), y{l - x" ~ y^)}, 

for example. 
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3. Expansions in Sobolev Orthogonal Polynomials 

Let H^{B'^) denote the space of functions for which (/, /)a is finite. This is 
not the space on i?** since the definition of (•,-)a require that / has second 
order derivatives. Nevertheless the standard Hilbert space theory shows that every 
/ e H^{B'^) can be expanded into a series in Sobolev orthogonal polynomials. In 
other words, 

oc cx: 

HHB'') = ^ ©V„^(A) : / = ^ proj„ /, 

n=0 n=0 

where proj„ : H^{B'^) i-^ V,^(A) is the projection operator, which can be written in 
terms of the orthonormal basis (|2.11|l as 

(3.1) proj„/(a;)= ^ HJ^ ^-^^lA^)^ flu = {f,Qlu)^. 

0<j<n/2 f=0 

where Hj = {Q^,^,QJ^„)a are independent of v as shown in H2.12(l . Let P^{x,y) 
denote the reproducing kernel of V^(A). In terms of the orthonormal basis 1)2.11(1 
in the previous section, the reproducing kernel can be written as 

0<j<n/2 f 

The projection operator can be written as an integral operator with P.^ as its 
kernel, which means that 



Proj„/(a;) = {f,P„ {x,-))a 
1 



A[{l-\\yr)f{y)]A[il-\\y\\)P^ix,y)]dy 



where A is applied on y variable. 

It turns out that the orthogonal expansion can be computed without involving 
derivatives of /. 

Proposition 3.1. For j > 1, let [3j ^ n - 2j + {d - 2)/2; then 

.oo^ fn _ 8j(j + 1) 

^ ' d^ vo\{B^) 

1 

"2 

furthermore, for j — 0, 



(/3,+j)(/3,+.7 + l) / f{x)Ql,{x)dx 



f{y')Yr''{y')du;{y') 



Proof By P;^^(Wo;x) = pj"'^^^ (2||a;||2 - l)Y^'-^^{x). Let j > 1. By Lemma 

Oand (ES, 

A [(1 - WxfQlA^)] = SjU + l)P]'jWo;x). 
Applying Green's identity 

, . AN, f f dv du \ , 

[uAv — vAujdx = / ~ 



dn d 



n 
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with v{x) = (1 — ||a;|p)/(.T) and u — Q"^ shows then 
8j(j + l) 



(3.3) 



Ad^ vol 
2j(j + 1) 



{l-\\x\\^)f{x)AP;^(Wo;x)dx 



d2 vol(B'') 



-2 



where we have used the fact that (1) — 1. Let dP^^'^^ denote the derivative 

of p(o./3). Using |(23 and it is easy to see that 

A[P;;^(iyo;a;)] =8 [2r2a2pjo.^.)(2r2 _ 1) 

+ (n - 2j + d/2)dPf'^'\2r^ - 1)] 

Let us denote the expression in the square bracket by Mj. The Jacobi polynomial 



(s) satisfies the differential equation 



(f - s^)y" - (-/3 + (/? + 2).s)y' + j(j + /3 + 1)2/ = 

Hence, changing variable 2r^ — 1 i-^ s, we conclude that 

1 
2 



2(1 - r^)M, - -j(j + + l)P,^"'''^'(s) + -(j + /3 + 1)(1 + 5)P,(lf +^^(s). 
On the other hand, using H2.8|l . (|2.7|l and the fact that 2, p. 782] 



(2j + /3 + 1)(1 + s)P,(lf +1^^^) = 2(j + [3)P'^'''{s) + 2jP,^if ^(s) 
we conclude that 

' "''^ ' ^ '-jpf'''^^(s) + (j- + i)/'|if^(s) 



j(l,/3)/ 



5(l:/3)/ 



2(l-.>/,..(^^+-^- + l)(^^+^-) 



2j+/3, + l 

= (/3, + J + l)(/3, + j)(l - r2)P,(!f ■)(2r2 - 1). 
Consequently, we have proved that 

(1 - \\xf)A[P-,{W,-,x)] = 4{(3,+j + l)(/3, + j)(l - r^)P^'J^\2r' - l)Yr''i^) 

= ^f3j + J + + j)Ql,ix) 

which leads to the stated result for j > 1 by (|3.3|l . The proof of j = is similar 
but easier, in which we need to use H2.13|l . □ 

Let us denote by P„(PF^; x, y) the reproducing kernel of V^(VFp), which can be 
written as 

P„(W^^;x,y)= ^ A~lP^{W^;x)Pa.{W^;y), 

\oL\—n 

where A^.p, = J ^j^[PaiyV^;y)YW^{y)dy in which is the normalization of W^. 
Let us also denote by C^(t) the Gegenbauer polynomial of degree n, and by a; • y 
the usual dot product of cc, y S R''. 
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Corollary 3.2. For f £ H'^{B'^) and x £ B'^ , 

proj„/(a;) = F„/(x) + (1 - \\xf) \ f f{y)Pn-2iW2;x,y){l - \\yr)dy 

(2) vol(B'*) Jfirf 

l<]<n/2 (1) 

where with x' = € S'^~^ , 

Proof. The values of Hj ~ {Q^i,, Q^i,)a are given in p.l2|l . It follows immediately 
that 

T.Ho'fo..QU^) = ^ I f{y')Y.Y:{y')Y::{x)d^{y')^YrJ{x), 
where the last step follows from the summation formula of spherical harmonics, 

v—\ i/=X ^ ^' 

for G 5*^"^. Furthermore, setting / — Q'j ^ with j >1 in (|3.2|1 also shows 



Hence, it follows from H3.2|) and l|2.12|l that 
(3.4) 

The relation 12.1411 readily shows that 

(3.5) JjQlMfdy - ^o\{B'')c2 jjPUW2;y)fW2{y)dy. 

We multiply (|3.4II by QJ^{x) and sum over and j. Using 1)3. 5|l and the fact that 
^j2,n~2j+ii^)^_^^ = j(j + l)/2, the stated result follows from (EU and JSH). □ 



It follows from this corollary that the orthogonal expansion of / with respect 
to (•, •)a coincides with the spherical harmonic expansion of / when restricted on 
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